Semiclassical Dynamics of Electron Wave Packet States with Phase Vortices 
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We consider semiclassical higher-order wave packet solutions of the Schrodinger equation with 
phase vortices. The vortex line is aligned with the propagation direction, and the wave packet 
carries a well-defined orbital angular momentum (OAM) hi {I is the vortex strength) along its main 
linear momentum. The probability current coils around momentum in such OAM states of electrons. 
In an electric field, these states evolve like massless particles with spin I. The magnetic-monopole 
Berry curvature appears in momentum space, which results in a spin-orbit-type interaction and a 
Berry/Magnus transverse force acting on the wave packet. This brings about the OAM Hall effect. 
In a magnetic field, there is a Zeeman interaction, which, can lead to more complicated dynamics. 
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Introduction. — The phase front singularities of a wave 
field, vortices, and related non-integrable phases have 
been introduced and examined in seminal papers [1-5]. 
Phase vortices appear naturally in the electron eigen- 
states in atoms, quantum Hall fluids, supermedia, ferro- 
magnets, Bose-Einstein condensates, and classical wave 
fields (e.g., in optics). While 2D vortices in condensed 
matter physics are point-like objects, with vorticity be- 
ing orthogonal to the plane of motion [6] , the optical vor- 
tices are mainly considered as linear objects in 3D space, 
with vorticity being aligned with the wave momentum. 
Wave beams with vortices constitute a fundamental set 
of modes with well-defined orbital angular momentum 
(OAM) [5]. Recently, these beams have found numerous 
applications both in classical and quantum optics [3,5]. 

Simultaneously, topological phenomena related to the 
semiclassical wave-packet dynamics of quantum particles 
have motivated intensive investigations in various areas 
of physics: condensed matter, high energy physics, op- 
tics, etc. [6-12]. This has resulted in the revisiting of 
the semiclassical equations of motion, but now with the 
Berry-phase terms taken into account, and the discovery 
of such phenomena as the spin Hall effect, with potential 
applications. In most cases, the topological Berry terms 
are due to the (pscudo)spin of particles, whereas the par- 
ticle's phase front is implicitly assumed to be locally of 
the plane-wave-type, i.e. without singularities. 

The dynamics of various types of vortex states with 
phase singularities also undergo the action of the topolog- 
ical Berry force, which can be associated with the Mag- 
nus force [13-15]. Both vortex and spin dynamics relate 
to such fundamental concepts as magnetic monopoles. 
Berry phase, space non-commutativity, and generalized 
Hamiltonian dynamics [1,2,4,6-16]. It has been shown 
very recently that the dynamics of optical beams with 
vortices reveals all the topological phenomena previously 
associated with spin, but now they are related to the 



intrinsic OAM carried by the vortex [15,17,18]. In con- 
trast to spin effects, these phenomena are polarization- 
independent and can be much larger in magnitude, since 
the OAM (the vortex strength) can take arbitrarily large 
values [15]. The OAM Hall effect (an analogue of the 
spin Hall effect) has also drawn attention in semiconduc- 
tor physics [19], but there it is also assumed that the wave 
packet states of the electrons have zero intrinsic OAM. 

While confined electron states with vortices are well- 
studied in molecules and 2D structures [20] , here we aim 
to analyze the semiclassical dynamics of a 3D propagat- 
ing electron wave packet or beam with intrinsic OAM 
due to a phase vortex. For simplicity, we will consider 
non-relativistic scalar electrons without spin. As we will 
show, electron wave packets with OAM manifest funda- 
mental topological and dynamical features which were 
associated so far mostly with spin dynamics. 

OAM states of a free electron. — The Schrodinger equa- 
tion in free space reads (^Hi-§i + 5^^^) — 0- Let 
us construct semiclassical (paraxial) wave packet solu- 
tions propagating along the z axis and characterized by 
a narrow distribution in (p, E) space around some cen- 
ter at (pc,i?c) = {pcez,pl/2m) . Usually, one assumes 
Gaussian-type wave packets with the maximal probabil- 
ity density and nearly-plane phase front in its center [8]. 
However, this is not the case for higher-order modes. By 
making the ansatz ip — exp ^ih~^ (p^z — -Ec^)] u, where 
u is a smooth function with respect to z and t, and ne- 
glecting the second-order derivative d'^u/dz^, we arrive 
at the parabolic-type equation 
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(1) 



Here 



+ Vl and d/dr = d/dt + (pc/m) d/dz 
is the time derivative in the coordinate frame {x,y,C,,T) 
moving with the wave packet center (C = z — pct/m, 
T ^ t). Since the operator in Eq. (1) is ^-independent. 
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such modes allow factorization, so that one can choose 
the transverse and longitudinal parts in the form of the 
known in optics Laguerre-Gaussian (LG) beams (with an 
azimuthally-symmetric intensity profile) and Hcrmite- 
Gaussian (HG) wave packets: 



/=-l 



.HG 



(C) 



(2) 



Here {r,f) are polar coordinates in the {x,y) plane, 
whereas / = 0, ±1,±2,... and m, n = 0,1,2,... are the 
quantum numbers corresponding to the azymuthal, ra- 
dial, and longitudinal directions. The zeroth mode, 
Z = m = n = 0, is a usual Gaussian wave packet. 

The explicit form of standard LG and HG solutions 
can be found in the optics literature [5], while here we 
will only emphasize their most significant features. First, 
they represent wave packets with Gaussian envelopes 
and constitute a complete orthonormal set of modes, 
so that (ui^m,n I"i',m%n') = ^ii'f^mm' ^nn' , and any local- 
ized wave function can be represented as a superposition 
\u) = ai^rn.n \ui.m,n)- Sccond, LG modcs with I ^ 

/,m,n 

contain a screw dislocation of the phase front on the wave 
packet axis, uf"^ oc exp{il(p); in other words, they con- 
tain a phase vortex of strength I at r = 0. Owing to this, 
the solutions (2) have a well-defined z-component of the 
0AM. Indeed, |wi,m,n) = hi |u;,m,n) {Lz = -ihd/dip), 
and, effectively, the electron possesses an intrinsic angu- 
lar momentum "L = h\ — Mez- The wave packets (2) 
also have a magnetic moment fi — gfiB^ [fJ-B — eh/2'm^ 
e = — |e|, c = 1), where 5 = 1 for classical orbital motion, 
but the factor can be different in general (e.g., g = 2 
for electron spin). 

The transverse distribution of the probability den- 
sity, p = |up, for LG modes with I ^ 0, represents 
(m -I- 1) concentric circles and vanishes at r = 0. At the 
same time, the probability current coils around z. j — 



'^[pPc + hlm{u*Vu)] 



^p{p^ + hle^/r), Fig. 1. 



(The same behavior is characteristic for the Pointing vec- 
tor of optical LG beams [5].) This implies that the elec- 
tron trajectories are effectively spiral in free space. This 
effect disappears in the classical limit h = 0, and can be 
regarded as a zitterbewegung due to the intrinsic 0AM. 
In what follows, we will consider only trajectories of the 
center of the wave packet (2), which is a sort of guiding 
center, propagating rectilinearly in free space. 

Dynamics in external fields. — In an exter- 
nal electromagnetic potential (A, $) generat- 
ing field (E,B), the Schrodinger equation reads 



e$ - J- (-ihV 

2m ^ 



ip — 0. The potentials 



„ , „,.v eA) 
are assumed to be smoothly space- and time-dependent 
(the fields are weak) to ensure the independent adi- 
abatic evolution of the modes (2). We assume that 
the wave-packet center, (pc,rc) (in what follows, the 
subscripts "c" are omited), moves in the vicinity of 
classical electron trajectory, p = eE -|- er x B, r = p/m. 
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FIG. 1: (Color online.) Transverse distribution of the proba- 
bility density p in LG beams with m = and different values 
of 0AM, /. Shown are the directions of common z component 
and different ip components of the probability current j. 



and approximately conserves its form (2) in the ac- 
companying coordinate system with the z axis locally 
directed along p {(^ = z — J pdt/m). This means that, 
effectively, we deal with a 'relativistic' configuration 
where 1 — Ip/p and the 'helicity', Ip/p = I, is conserved 
(see below). The local coordinate frame is transported 
along a curved trajectory and. This produces a phase 
shift due to the Berry phase [4] and a deflection of 
the wave packet center due to the geometrical force 
[6-12,15], both described by the effective Berry gauge 
potential (connection) and field (curvature). 

Geometrically, the Berry connection provides for the 
parallel transport of the state vector over the phase 
space [4,7]. Since the z axis is attached now to p, the 
solution (2) becomes essentially momentum-dependent. 
The parallel transport of the wave packet over momen- 
tum space implies the covariant derivative D/Dp = 
d/dp + A'-'^p), where A^'^ = i {ui,„,,,\d/dp\ui,^,„) is 
the Berry connection, whereas the corresponding curva- 
ture is B*-'^ = d/dp X ^'■'^ . As will be seen, only the quan- 
tum index / noticeably contributes to the Berry connec- 
tion. Furthermore, {ui^m,r\\d/dp\uir^rn',n' mm' f^nn' ; 
which evidences the independent evolution of the modes 
(2) . As in the relativistic case of massless particles with 
spin, the Berry gauge field takes the form of a 'mag- 
netic monopole', which originates from the local vor- 
tex structure exp(ilip) in the wave packet [15]. Indeed, 
exp{il(p) — + ieyf , where e is a unit vector orthog- 
onal to p. Under variations of p, e moves on the unit 
sphere p/p, which leads to the magnetic-monopole-type 
connection A. ^ i {cx — iey) d/dp {cx + iCy) and corre- 
sponding curvature B = —p/p^ [4,11,12]. As a result 
we have ^^'^ = lA and B^'^ = IB, so that each mode is 
characterized by the 'charge' I in the 'magnetic monopole' 
field in the momentum space. 

The semiclassical dynamics of the wave-packet center 
can be described involving the minimal coupling pre- 
scription for the electromagnetic and Berry's gauge fields. 
This results in the Lagrangian [8] 



C= -'H + pr + eAr + hlAp . 



(3) 



Here the Hamiltonian acquires the energy correction due 
to the Zeeman interaction of the intrinsic 0AM with the 
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magnetic field [8,9,12]: n^^ + e<^> + A = E, 
— /xB, where the magnetic moment fi = gfJ-B^ is also 
momentum-dependent now. The Berry phase term, last 
in Eq. (3), is of the form of the relativistic spin-orbit 
interaction [12], but now this is the orbit-orbit interaction 
between the intrinsic 0AM and the external degrees of 
freedom. Eq. (3) yields the common phase of the wave 
packet : 

e = fr^ j (pdr - Edt) + tr'^e j Adr + I J Adp , (4) 

which substitutes for the free-space one, {pz — Et). 
The three terms in Eq. (4) are, respectively, the dynam- 
ical phase, the Dirac (Aharonov-Bohm) phase [1,2], and 
the Berry phase [4,7,8]. The latter provides for the par- 
allel transport of the transverse structure of the wave 
packet along the curved trajectory [15,17]. 

Considering the Euler-Lagrange equations for C — 
£(p, p, r, r), Eq. (3), we arrive at the semiclassical equa- 
tions of motion for the wave-packet center [6-16]: 

^ /\ 

p = eE - — + er X B , f = — + hip x B . (5) 

or m op 

These equations follow from the Hamiltonian formalism 
as well, where the minimal coupling with the electro- 
magnetic and Berry's fields implies a deformation of the 
symplectic structure ft — g^^dXi A dXj/2 [7,16]: 

ri = dpi A dri + e-^Bkdri A drj + M-^Bkdpi A dpj . 

Here X = (r,p), and g"^^ is the symplectic metric ten- 
sor. The equations of motion (5) take the simple form 
g^-'Xj = dH/dXi, whereas the Poisson brackets of the 
dynamical variables (or commutators of the correspond- 
ing operators) become non-trivial, {Xi^Xj} = g^ (cf. 
[7,12,16,21]): 

{Pi,P]} = D^^efijkBk , {ri,rj} = D^^hhijkBk , 
{r„pj} = D-^ {6,j - ehlB.Bj) , D^^deig^^ .(7) 

Here D = \ — ehlUB is a correction to the phase space 
volume which modifies the density of states [21]. 

Let us briefly analyze the most remarkable features of 
equations (5). The last term in the second equation (5), 
given by hip x p/p^, represents the anomalous velocity 
(f i/( p) or the 'Lorentz force' from the Berry's 'mag- 
netic monopole' in momentum space [6-12]. It causes the 
transverse deflection of the wave packet trajectory, which 
is proportional to the intrinsic 0AM hi [15]. This results 
in the transverse orbital current, or 0AM Hall effect, 
similar to the spin current in the spin-Hall effect [10-12]. 
For instance, when B = 0, E = const, Eqs. (5) take the 
simple form: p — eE, r — p/m -f- ehVE x p/p^ , and can 
be readily integrated analytically (see [10]). They show 
the transverse shift of the trajectories, equaling hl/pQ 



when po-LE, Fig. 2. The 0AM Hall effect can be much 
stronger than the spin one, since it is proportional to 
which formally can take arbitrarily large values. At the 
same time, the anomalous velocity is a counterpart of 
the Magnus force for 2D vortices which also appears due 
to the Berry curvature term, but in coordinate rather 
than momentum space [13,14]. In an external magnetic 
field, E = 0, B = const (for simphcity, let B ± p), 
Eqs. (5) are reduced, in the linear approximation in h, 
to p = ^p X B, r = m"^[p + ehl{l ~ f )B/p]. Thus, 
at 5 7^ 2, an OAM-dependent transport of electrons can 
appear in a constant magnetic field. 




FIG. 2: (Color online.) Central trajectories, Eqs. 5, of the 
electron states with different values of 0AM, Z, moving in a 
uniform electric field E = E&y (po = poCz, ro = 0, and / = 
line corresponds to the classical trajectory). 

The intrinsic 0AM, 1, can also be considered as an in- 
dependent intrinsic dynamical variable. It can be shown 
that the helicity Ip/p is conserved during the evolution, 
Eqs. (5), if 1 obeys the equation 

i=-(-B + ^)xl. (8) 

\m P ) 

Eq. (8) is of the form of the well-known BMT equation 
for the electron spin precession [22,12], where the Berry 
phase (spin-orbit) term is taken in the relativistic limit 
(because intrinsic 0AM reveals topological features sim- 
ilar to the spin of a massless particle), while all other 
terms are non-relativistic. Eq. (8) is compatible with 
Lagrangian (3) only if either g = 2 or B = 0. Under 
evolution of the 0AM with g 7^ 2 in a magnetic field, 1 
does not follow p (the precession frequency of 0AM dif- 
fers from the cyclotron frequency), and thus the initial 
assumption of this study is not valid. In this case, the 
problem becomes much more complicated and requires 
an independent study. 

Discussion. Comparison with optics. — Electron states 
with vortices can appear in both bulk solids and in free 
space. These may open new avenues of research in con- 
densed matter physics (0AM Hall effect) as well as in 
electron microscopy and holography [23] (similarly to sin- 
gular optics). A comparative analysis of the electron 
0AM states and their optical counterparts is summarized 
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in the Table I. Dynamics of optical beams with vortices 
is described by equations similar to Eqs. (5), where the 
refractive index of the medium, n, plays the role of the 
external scalar potential <I> [15]. At the same time, there 
is no counterpart of the external curl potential and mag- 
netic field in optics. Therefore, all electric-field-related 
features arc absolutely similar in electron and optical sys- 
tems, while the magnetic-field-related effects are inherent 
to the evolution of charged particles only [24]. 

It is important to find some ways for creating elec- 
tron phase vortices. On the one hand, the phase vortex 
cannot be created from a plane wave by a large-scale 
electromagnetic potential. Indeed, the phase incursion 
at the closed contour around the vortex line equals 27r/. 



A Dirac phase of this kind can appear only due to a real 
magnetic monopole [1,23] or in the presence of an infi- 
nite Aharonov-Bohm solenoid [2]. On the other hand, 
a phase vortex can be created by a short-scale poten- 
tial. For instance, an optical vortex can be produced us- 
ing a diffracting grating with an edge dislocation ('fork'). 
For electrons, such a grating can be provided by a real 
thin crystal plate with a dislocation. This offers the op- 
portunity of using phase vortices in electron or neutron 
crystallography. One can also use a bulk solid as an ef- 
fective refractive medium for free electrons, because the 
electron changes its effective mass there [23]. Then, a 
spiral-thickness solid plate will create an electron phase 
vortex similarly to the analogous optical lens [3,5]. 





Optics 


Electrons 


wave field 


E 




external potential field 


n, Vn 


E 


external curl field 




A, B 


generation of 


hologram (grating with dislocation) 


crystal plate with dislocation 


phase vortex 


spiral-thickness lens 


spiral-thickness plate 






magnetic monopole 


common features 


orbit-orbit interaction: Berry phase and Magnus/Berry force 


distinctive features 


Lorentz force, Dirac phase, Zeeman interaction [24], modified density of states 



Table I: Comparison between electron 0AM states and their optical counterparts. 
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